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Abstract. In quasicrystals, there are not only conventional, but also phason displacement fields and asso- 
ciated Burgers vectors. We have calculated approximate solutions for the elastic fields induced by two-, 
three- and fivefold straight screw- and edge-dislocations in infinite icosahedral quasicrystals by means of 
a generalized perturbation method. Starting from the solution for elastic isotropy in phonon and phason 
spaces, corrections of higher order reflect the two-, three- and fivefold symmetry of the elastic fields sur- 
rounding screw dislocations. The fields of special edge dislocations display characteristic symmetries also, 
which can be seen from the contributions of all orders. 

PACS. 61. 44. Br Quasicrystals - 61.72.Lk Linear defects: dislocations, disclinations - 62.20.Dc Elasticity, 
elastic constants 



1 Introduction 

It has been demonstrated frequently (see, e.g., j^]) that in 
quasicrystals, as in periodic crystals, plasticity is caused 

• by dislocations moving under external stresses. If the elas- 
' tic fields around dislocations are known, their influence on 
. the lattice and the interaction between dislocations can be 

■ calculated. In quasicrystals, dislocations are surrounded 
not only by phononic, but also by phasonic elastic fields. 
They are described by a generalization of the standard 

' elastic equations. The generalized elastic theory has be- 
come a powerful and important tool for studying the me- 
chanical behaviour of quasicrystals. 

The dislocation problem is solved when the medium's 

• elastic Green's function is available. The elastic Green's 

■ function has been established in closed form for isotropic 
, 1^ and hexagonal |^ ordinary media, and for pentago- 
nal 01, decagonal and dodecagonal § quasicrystals. An 
approximate solution for the elastic Green's function of 
icosahedral quasicrystals is given in In this paper, the 
results are obtained by direct solution of the equations of 
balance instead by use of Green's method. 

Analytical solutions for the elastic fields around dis- 
locations exist only for the above mentioned cases 
1^. In this paper, we present the solution for icosahedral 
quasicrystals in terms of appropriate perturbation series. 

The structure of the paper is as following. Starting 
from the density wave picture, we summarize in Section 
^ the generalized elastic theory of icosahedral quasicrys- 
tals, including dislocation elastic theory. Section ^ deals 
with the generalized projection method as one possibility 
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to solve the dislocation problem. Here we present a set of 
recursion formulae, which can be used to calculate pertur- 
bation expansions of the elastic fields. In the last section, 
we discuss displacement fields induced by different types 
of fivefold dislocations in icosahedral quasicrystals. 

Most concepts and notations in this paper are identical 
to those used in Q . 

2 Elastic theory of icosahedral quasicrystals 

2.1 Some fundamentals 

A quasicrystal is a translationally ordered structure with 
sharp diffraction pattern exhibiting non-crystallographic 
symmetry. For this reason, its mass density p(x) can be 
written as a sum over density waves: 

p(x) = P^^"^'^ - E IPk|e*('^-''+*'') . (1) 

keL keL 

Here, L is a module over the reciprocal quasilattice. The 
numbers (/)k are the phases of the complex coefficients pk- 
For icosahedral quasicrystals, the diffraction pattern 
and L display icosahedral point symmetry. Six vectors k^j , 
a — 1, . . . ,6, pointing to appropriate six of the vertices of 
an icosahedron can serve as a basis of L. Phenomenolog- 
ical Landau theory shows that there are six degrees 
of freedom, which can be interpreted as the phases 
of the basis vectors. Thus, a frequently used approach is 
the extension of the density (|^) to the density of a peri- 
odic structure in six-dimensional hyperspace, which can 
be subjected to a six-dimensional displacement field 7. 
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is recovered by a cut of the hyperperiodic structure by 
physical space Q . A parametrization of the (^k is 



Ok,o 



K • 7, 



which involves vectors k belonging to the reciprocal hy- 
perlattice in six dimensions. The set {ho} must be linearly 
independent over the real numbers and serves as basis of 
the reciprocal hyperlattice, which is usually chosen to be 
hypercubic. Then, the direct hyperlattice is hypercubic 
also, with basis vectors /?=!,. ..,6, defined by 



(3) 



The icosahedral group Y acts on the hyperspace ac- 
cording to the reducible six-dimensional representation 

= . The hyperspace decomposes into two or- 

thogonal, three-dimensional invariant subspaces i?" and 

, belonging to the irreducible representations (vec- 
tor representation) and of Y. i?" is the physical or 
parallel space and E-^ the perpendicular space. Two pro- 
jection operators P" and P-*- can be applied to hyperspace 
vectors to obtain their respective components in i?" and 
E^. This procedure is shown in Figs. ^, || for the ortho- 
normal natural basis {e^j}, spanning a hypercube in six- 
dimensional hyperspace. 

All quantities depend on the physical space coordi- 
nates x" — X only. Because of the orthogonality of ij" 
and E-^ , Eq. (0) can be written 



= 0k, - kll • u - k-^ 



(4) 



Notations k = k" ©k-*- and 7 = 7" ©7"'" = u©w are used, 
where k" = k £ L. u is the ordinary phonon displacement, 
whose character is propagating, w is the phasonic dis- 
placement with diffusive character | ]lO| . In the atomic pic- 
ture, non-vanishing w-fields lead to local rearrangements 
of atoms, which are called phasonic flips. 



Spatially varying 7 and 0k, respectively, belong to de- 
formed states, which can be described by the elastic tensor 
fields of distortion (3 and strain e defined by 



du = /3"dx , 
.. 1 



dw = /3™dx , /3 = 



/3" 

£"" 



(5) 
(6) 



In case of single- valued displacement fields u and w. 



duj 

1 / dui duj 

2 \dxi dxi 



r = ^ (7) 



where i,j G {1,2,3}. The six symmetric and three anti- 
symmetric components of /3™ mix under the action of Y, 
and therefore no invariants can be constructed from the 
symmetric or antisymmetric components of (3^ only. So 
due to (^, the phason strain tensor e"-' must be the full 
/3™ to contribute to the elastic energy density F. 

The relations (I) and d) hold in the linear regime 
ifx^l' ifi^l ^' "^tiich we are interested in. Linear elas- 
ticity is described by an elastic energy density which is an 
Y-invariant quadratic form of the components of e: 



(9) 



Here, the Greek indices a,(3 £ {I,.-. ,6} refer to 
(a,/3 £ {1,2,3}) and to E-^ (a,/3 £ {4^,6}). This is 
different from the meaning of a, /3 in Eq. (|3|). 

Differentiation of leads to the generalized Hooke's 
law: 



dF 

dr ■ 

L/C 



C£. 



(10) 
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with the generahzed stress field 



a" 



(11) 



where <t" — cr'"'* is symmetric. The stress tensors cr"-'", 
applied to normal vectors n of real or fictitious surfaces 
in physical space, lead to surface forces t"'"" = (T"'™n, 
which must be applied to keep the system in balance. 
These forces have components in i?" and E-^, respectively, 
and are combined to t t" t'" Besides the con- 
ventional body force f", an analogous force f™ with direc- 
tion in E-^ has to be introduced to formulate the general- 
ized equations of balance consistently. We use the notation 
f = p©f"'. 

Hooke's law can be formulated very simply with the 
help of group theory. Due to Eqs. (p|), e" transforms as 
the representation {r^^ E^)sym = 1^(BE^ under actions 
of Y, whereas e'" transforms as E^ (g) E^' = 7"^ ® T^. The 
transformations of the components of cr"^™ are quite the 
same. With the irreducible strain components related to 
the coordinate systems of Figs. |l| and ^ 



U _ \ U It 

^1 - y^ri\ 

^■w _ \ w 



[el 



r5, . 



W It w _ \ w 
, t.pi\ ,£.5 — [t p5 , . 



It 



(12) 



They can be measured by ultrasound transmission [ p^ . 
The elastic constants /14 and /i5 belong to pure phason 
elasticity. They are determined currently onl y in directly 
from diffuse scattering around Bragg peaks Q. Values 
for the phonon-phason-coupling /X3 have been calculated 
in computer simulations to have much smaller absolute 
values than the two phononic elastic constants [p^ . 

Isotropic phonon elasticity in thermal equilibrium re- 
quires decoupled phonon and phason elasticity, i.e. /Z3 = 0. 
Isotropy in phason elasticity is given in the spherical ap- 
proximation /i4 = /Z5 discussed in in addition to the 
condition /i3 = 0. 



Gauss' theorem, applied to Eq. (10), provides the well 



known elastic equations of balance in generalized form: 



div (T + f = , 



dxi 



(16) 



where /i = . . . , /g = With the help of Hooke's law 
and (||) , div cr can be written in terms of the displacement 
field 7: 



div(T = D(V)7, (divcr)„ = D„^(V)7/3 , 



(17) 



(see 
ing t: 



12 and Appendix A) and analogous vectors contain- 
le irreducible stresses, Hooke's law reads 



>i 
/12 ^3 
//4 

. /13 ^5 



£1 

^5 
£4 



(13) 



Since only equal-indexed components of the same irre- 
ducible representation can interact with each other, the 
number of independent second-order elastic constants is 
restricted to five. The elastic energy density correspond- 
ing to (O) is 



E 



2 ^4 £4 



2 £5 



£5 



-5 • 



(14) 



A comparison with (^ yields the coefficients Caifij- In 
Appendix C, we discuss the conventions for independent 
elastic constants used by other authors. 

Mechanical stability requires F > for every £ ^ 0, 
which is fulfilled when all eigenvalues of the elastic tensor 
of ( p^ are positive. This leads to the following constraints 
on the elastic constants: fii > 0, /i2 > 0, /i4 > 0, /i5 > 0, 

and pi2M5 > Ms- _ _ 

According to ( |l3|) and (|l4|), the elastic constants /ii 
and /i2 describe pure phonon elasticity. Without phason 
elasticity, for example on a short time scale, icosahedral 
quasicrystals behave like isotropic media with the two 
Lame-constants 



dxidxj 



(18) 



Here, 71 = Ux, ■ ■ ■ , 76 = Wz- According to (|17|), the equa- 
tions of balance (Qq) result in 



D(V)7 + f = 0. 



(19) 



Eq. (|19[), f = PQf^ and 7 = u©w imply a decomposition 
of D(V) into four 3x3 blocks: 



D(V) 



Explicitly, we haveQ 



D™(V) D™(V) 
D'""(V) D'"'"(V) 



D™(V) = 1 -f (A -I- Ai)V ® V , 
D'™(V) = D"'"^*(V) 

Ei{x,y,z) E3{x,y) E2{x,z) 



M3 

V6 



F2{y,x) Ei{y,z,x) Fz{y,z) 
E3{z,x) E2{z,y) Ei{z,x,y) 



D™"(V) =Ai5l V 

, M4 - A*5 



Ei{x,y,z) E^{x,y) E^{x,z) 
Fb{v,x) Ei{y,z,x) E;i{y,z) 
F5iz,x) Ez{z,y) Ei{z,x,y) 



(20) 



(21) 



^ = ^(Mi - M2) . 



(15) 



^ Here, some printing errors of Ref. ||q| have been eliminated. 
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In (21), we have used the abbreviations 



Fi(a,6, c) 
F2(a,6) 



92 1 92 
2 52 



(9c2 



F3(a,6) = 2t 



92 



dadb 



(22) 



F4(a,5, c) 



a&2 



dadb 



T = i(l + -\/5) is the Golden Mean, defined as the positive 
root of the quadratic equation — x — 1 = 0. 

2.2 Green's function of the elastic equations of balance 

The solution of (19) for arbitrary body forces f can be 
obtained by calculating the integral 



7(x) = G(x-x')f(x')d3x', 



(23) 



where G(x — x') is the elastic Green's function of icosa- 
hedral quasicrystals. The constitutive equation for G(x) 
is 



D(V) G(x) + 1 (5(x) = . 



(24) 



An approximate solution for G(x) is given in Since 
(21]) are the components of D(V) for linear elasticity, the 
displacement field provided by ( p3| ) is the correct solution 
only in domains where the linearity conditions are fulfilled. 

For the purpose of this work, we only need to know the 
solution Goo (x) for elastic isotropy in phonon and phason 
spaces, defined by the conditions ^13 — and /i4 — /i5. 
With the well-known elastic Green's function for three- 
dimensional, isotropic media and the solution of the 
fundamental Poisson's equation in three dimensions, the 
exact Goo(x), which reflects decoupled phonon and pha- 
son elasticity, is 



G?f?(x) 
G-(x) 

GoT(x) 



1 



87r^(A -t- 2^) 
1 



{X + 3fi)-^l + {X + 



47r^5|x| 



Go"ir(x) 



(25) 



2.3 Elastic theory of dislocations 



A dislocation D is characterized by its non- vanishing line 
integral 



Ak = 27rmk, mk = 0,±l,±2. 



'dFo 



along any closed contour OFd surrounding the core of D, 
which exists in physical space only [^|l^. Eq. (|6|) guar- 
antees for the continuity of the mass density /^(x) outside 
the dislocation core, as can be seen from (|^). In case of 
periodic crystals, the lattice remains unaltered outside the 
dislocation core. Because of phasons, the same is not true 
in quasicrystals . 

From Eq. (||) , condition ( p6|) can be expressed in terms 
of K and d7: 



K, ■ d-y — — 27r mt . 



(27) 



The six-dimensional Burgers vector is b = b" b™ , where 
b", b™ and b are defined 



b" = * du , b"' * dw , b = * d7 . (28) 

JdFo JdFn J dFo 



This allows us to rewrite (^^ in the form 
K • b = k" • b" + k-L • b"" = -27r 



(29) 



As a conclusion from (^ and (^9|), b must be a vector of 
the direct hyperlattice. The irrational orientations of iS" 
and E-^ in hyperspace provide non- vanishing components 
b" and b™ for every dislocation. 

Because of Eqs. (p8|), the displacement fields of dislo- 
cations are multiple-valued, and therefore not integrable. 
According to potential theory, the non-integrability of Eqs. 
(||) corresponds to curl /J"'"" 7^ 0. From Stokes' theorem, 



d(u, w) 



OFd 



/3"^'" dx 



dFu 



curl/3"''"df = / a"'"'df, (30) 

Fd JFo 



where the surface integrals extend over the surface Fu, 
bounded by dFn and pierced by the dislocation D. Be- 
cause of Eq. (pO[), the tensors a"'*" are dislocation den- 
sities in phonon and phason spaces per area of physical 
space. It is 



curl/3'' 



dxk 



(31) 



Eqs. (|31|), in connection with boundary conditions to be 
fulfilled, are the defining equations for /3"'™. 

By means of Eqs. (^, the elastic distortion fields /3"'™ 
lead to the strain fields e"'™. Connected with the elastic 
strain fields of dislocations are stress fields ( [lO|). In the 
absence of body forces f , they must fulfill Eq. (|16[) in the 
form 



div (7 = 0. 



(32) 



3 Solving the problem of a single dislocation 



(26) Elastic fields fulfilling Eqs. (|l^, (H) and (||) simultane- 
ously constitute the solution of the dislocation problem. 
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Note that Eqs. ( |28| ) for a single dislocation are a special 
case of the general situation, which is characterized by 
an incompatibility field The projection method |Q 
is a possibility to solve the dislocation problem. A sum- 
mary, containing the description of several other methods, 
is given in 



3.1 The projection method 

The idea of the projection method is the following. In a 
first step, an arbitrary, multiple- valued displacement field 
7 = u © w fulfilling (Gq) , 



du , 



9Fc 



dw . 



d7, (33) 



dFr- 



has to be found. According to Eq. (H), this displacement 
field leads to distortions, and from B) the corresponding 
eigenstrains e"'™ can be calculated. But the eigenstress 
field (T = C e will not be compatible with the equations 
of balance (|3|) in general. There may exist a divergence 
divo- ^ 0. 

In the second step, a single- valued displacement field 
7 = u ® w, leading to strain fields e"'™ and stresses & = 
C e, due to Eqs. (||) and (|o|), respectively, has to be found. 
7 has to be determined in such a way, that the divergence 
of the total stress field cr — & + & vanishes: 



div (T = div {& + (t) 



0. 



(34) 



In this step, the stress field & is projected from the space 
of stress fields with arbitrary divergence onto the space 
of stresses with vanishing divergence. The result of this 
projection is the true stress field cr. Since the displacement 
field 7 is single-valued, we have 



du = , * dw = , * d7 = . (35) 

dFn JdFo JdFo 

Finally, the elastic fields of the dislocation are 
7 = 7 + 7j £ = £ + £, (T^a + a-. (36) 

Eq. (|l0| ) is satisfied because we have used {a, &) = C (e, s). 
(||) is fulfilled because of Eqs. (|3|) and (||). (|2|) is true 
due to Eq. (p^. 7 is multiple- valued because 7 is, whereas 
£ and <T are single- valued. 

The problem of finding the right 7 can be solved as 
described below. The equations of balance (|9|) read for 
the case of the displacement field 7 



D(V)7-Kf = 0, 

where D(V)7 = divd-. Comparing 
leads to the identification 

f = div IT 



for the right f to be used as fictitious body force in ( p7D 
in order to calculate 7. Therefore, in the second step of 
the projection method, the equations of balance 



(37) 

with Eq. (|4|) 
(38) 



D(V)7 + divo- = 



(39) 



have to be solved. We want to stress that the body force 
f — div (T is not a real force. 

The projection method amounts to adding a singular 
displacement field and a non-singular one in order to sat- 
isfy Eq. (p2|), which is not yet fulfilled after the first step. 
The splitting up into a singular and a non-singular part 
and the introduction of local 3jXGS ; clS done below, has al- 
ready been discussed in |l6[. 



3.2 Symmetry-adapted coordinate systems 

In this paper, we consider straight dislocations D having 
two-, three- and fivefold line directions. For simplicity, a 
coordinate system with its z^-axis parallel to the re- 
spective dislocation line of D should be used in physical 
space. There arise certain symmetries from phasonic Burg- 
ers vectors, which can be seen clearly when additionally 
choosing the ^"'-axis of the new coordinate system in 
phason space parallel to that symmetry axis, which corre- 
sponds to the line direction in physical space. 

The difference between the two-, three- and fivefold 
dislocation lines lies in different components of the elas- 
tic tensor C. The components Caipj resulting from (H), 
( |l^ belong to the coordinate systems of Figs. |l|, ||, which 
are appropriate coordinate systems for twofold disloca- 
tions. The components CD,aii3j of the elastic tensors in 
symmetry-adapted coordinate systems for the three- and 
fivefold dislocations arc obtained by coordinate transfor- 
mations. 

Consider orthogonal transformations in phonon and 
phason spaces. 



(40) 



connecting the components of the vector v in the symme- 
try-adapted coordinate systems with its components be- 
longing to the coordinate systems of Figs. ^. The trans- 
are then 



formations of the strain components are 



^D,ij — ^ D,ik^ D,jl^kl J 

Explicitly, we have taken 



'D,ij 



— ^ D,ik-^ DJl^kl 



(41) 







I 









I 


3 


= Cl 




-I 


' 3 


= Cl 


I 






I 








1 






"t - 


-l" 






" I Or' 


5 


= C2 


I 





' 5 




I 






I 


r 






-r I 



(42) 



as orthogonal transformation matrices in ( ^ ) . They result 
from choosing the directions P"(e3 + 64 — eg) and P"e3, 
respectively, as dislocation lines in Fig. |l|. The coefficients 
Cl, C2 are 



1 



Cl 



C2 



(43) 



When substituting the components e^"- of Eq. (^) by com- 
ponents e^jT/' '^^'^ ll^*^ inversion of Eqs. (^), the com- 
ponents CD,aif5j are obtained. 
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The positions in phonon space are best described by 
cylindrical coordinates r and (f), beside the cartesian coor- 
dinates X, y in the symmetry-adapted coordinate systems. 
From now on, we use the notation x = [x, y]* for the po- 
sition perpendicular to the dislocation lines. 



3.3 Generalized differentiation and its application to 
the projection method 

The best 7 to start with in the first step of the projection 
method is 



The resulting distortion and strain fields are 



(44) 



= /3 ' = 



1 



y 



^2 y h' 2; 

y h 2; 



47r(a;2 



-2b^y b'lx-b^y -h^y 
b\x - b^y 2 b^x 



b^x 



blx 




(45) 



The corresponding dislocation densities a"'™ ( pl| ) are 
for X ^ 0. To consider the behaviour at x = 0, where 
no partial derivatives in classical sense are defined, one 
has to assume all functions to be generalized functions 
and calculate the generalized partial derivatives 1 22 . Note 
that, because of (\ 



and (^), the linear elasticity breaks 
down below a finite cut-off distance from the dislocation 
line. 

In case of homogeneous functions / of degree —1 in 
two variables x, y, the generalized partial derivatives are 





dl 


dx 


dx 


9/ ^ 


dl 


dy 


dy 



+ 5{^)i fdy, 

G JOG 

~ f fdx. 

G JdG 



(46) 



where the respective first expressions represent certain 
generalized functions connected with the classical partial 
derivatives. G is any bounded region containing the origin 
X = 0. Eqs. (^6|) indicate that there may be (5-functions 
on the dislocation lines. Their respective weight may de- 
pend on the choice of G, but becomes independent of G 
when considering any combination of generalized partial 
derivatives with the total classical derivative vanishing. 
Calculating the dislocation densities ( ^l|) results in 



^ ^^^^ 



j-iuy^j y lAi' -1-/3,2' 

dG 

u.w 



= 5 



3j3 



2tt 

' U.IU 



-^6(^)6 d0=5,3&r'5(x) 

JdG 



Therefore, the dislocation densities are point-like, 
a"^"" ^ (5(x) [0, 0, b"^"-] . 



(47) 



(48) 



According to (M), a"'"" = 0. 

From Eq. ( [lO| ) and Eq. (^5|) , the stress field & entering 
the projection method in the second step is homogeneous 
of degree —1. Therefore, Eqs. (^) apply when calculating 
div &. 

In the second step of the projection method, the solu- 
tion of ( ^ ) is split into further two steps. In a first step, 
( ^9|) is solved for x 7^ 0. This means that only the classical 
derivatives {div &)\:x^o enter (^9|). Denoting the solution 
7', Eq. (p9) becomes for x 7^ 



D(V)7' -hdivo- = 0. 



(49) 



The displacement field 7' leads to the strain field e' and 
the stress field cr', due to ^ and (p^). The divergence 
of the stress field & + cr' will vanish when considering 
classical derivatives only, but may still have components 
proportional to 5{'x.). In a second step, these point-like 
stress sources have to be compensated. With div {a+cr') = 
(5(x)f£', Eq. (M) becomes 



D(V)7" + '5(x)fi' = 0. 



(50) 



Here, f£' is a constant vector having the unit of a force per 
length. The solution 7" leads to the strain field e" and 
the stress field cr" . The final results entering Eqs. (|36| ) are 



7 = 7 



7 



(51) 



In a short excursion, we want to discuss the meaning 
of the divergence (5(x) f£' of the stress tensor on the dislo- 
cation lines. According to (p6[), 

/La=/ [(^ + ^')aldy-(^ + <T')a2dx] . (52) 
JdG 



(p2[) is the a-component of 



f t" ds = f (a- + cr')nds , 
JdG JdG 



(53) 



where t" is the six-dimensional surface force remaining af- 
ter the divergence of the stress tensor in terms of classical 
derivatives has been brought to zero, t'l represents the net 
force per length of the z"-axis acting on the interior of G, 
which must vanish [E3|. 



3.4 Recursion formulae 

Unfortunately, there is no possibility to calculate 7 in the 
second step of the projection method exactly. We have ex- 
panded the elastic fields into perturbation series provid- 
ing an approximate solution of the dislocation problem, 
the zeroth order being the solution for the case of elastic 
isotropy and the higher orders the perturbation arising 
from the deviation from isotropy. Closely following the 
projection method, the solution of the recursion formulae 
is the summation of the singular displacement field of ze- 
roth order and non-singular displacement fields of higher 
order. This kind of dealing with anisotropy has already 
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been used in connection with Green's functions of crys- 
tals and quasicrystals Q and to calculate the fields 
of dislocations and cracks |25). ||2^ comprises a more 
enclosing approach to the handling with elastic anisotropy. 

The zeroth order is the solution for elastic isotropy 
(/i3 = and /i4 — /is), which is known (see Appendix 
B). The elastic fields can be expanded with respect to 
the variables fj,3 and — fx^. We shall mark any physical 
quantity / being proportional to fi^i^i^ — ^5)" by lower 
indices mi^rnn Ai"(Ai4 - Ms)"- 

From (pl|), it is obvious that 



D(V) = Doo(V) + Dio(V) + Doi(V) , 
where according to Eqs. ( pT| ) 



Doo(V) = 



/.i 1 + (A + /i)V «) V 

/islV^ 



(54) 



(55) 



Because of (p^, the decomposition ( |54| ) corresponds to 
the following decomposition of the clastic tensor C: 

C - Coo + Cio + Coi . (56) 
For the elastic fields, the perturbation series read 

The conditions, which the different orders must fulfill, 
are the following. First, all displacement fields except for 
zeroth order must be single- valued: 



h = f d7oo : 

JdFo 







dFc 



d7„„ (mn^OO). (58) 



This is because b belongs to zeroth order. Second, Hooke's 
law (no) must connect strains and stresses in all orders: 



""00 = Coo £00 I 

CmO — Coo £mO + Cio £m-l,Q , 
CfOn = Coo £Ori + Coi £o,ri-l , 
mn — Coo £mn + Cio £m-l,n + Coi £m,n-l ■ 



(59) 



Here, the last three equations are valid for m > 0, n > 
and both m, n > 0, respectively. Third, the divergence of 
the stress field of each order must vanish: 



(60) 



The constitutive equations (58), 1)^^ and (|60|) can be 
solved recursively using the projection method. At this, 
the zeroth order must be treated somewhat different than 
higher orders. 

Note that 7 and £ given by Eqs. (^4|) and ( ^5| ) belong to 
zeroth order, as the stress field &qo = CqqS. divcroo ~ 
must be valid, and therefore the projection method re- 
quires the solution of the equations of balance 



Doo(V)7oo +div(Coo£) = 



(61) 



When omitting th e lo wer indices, this is exactly the pro- 
cedure of Section 3.1. The solution of (^ij) leads to the 
strain field eoo according to Eqs. (||). Finally, we have 

7oo — 7oo + 1 T ^00 — £00 + £ , Coo — Coo (£00 + £) ■ 

(62) 

The full 7oo is given in Appendix B. 

We demonstrate the procedure in case of higher orders 
for toO 7^ 00 only. The approach in the other cases is 
completely analogous. When calculating the elastic fields 
of order mO, the solutions of all orders rhO with rh < m 
have already been determined. The condition for the stress 
field (Tmo to be fulfilled is divcroo = 0, where a^o is 
determined by the second Eq. (^). Obviously, the solution 
of the equations of balance 



Doo(V) 7„o + div (Cio £m-i,o) = 



(63) 



provides the displacement field 7^0 of order mO. Eq. ( |6^ ) 
is obtained immediately by the action of the differential 
operator div on both sides of the second Eq. (|5^) and the 
fact that Doo(V) 7^ = div (Coo £mo), with and £,nO 
connected by Eqs. (g). 

Concerning the perturbation exp ansions of the fields 7, 
e, & and f introduced in Section 3T , we have 7„u, = 7m„ 
and £mn = £mn except for zeroth order, where (|62|) must 



be applied. The relation cr^ 



is true for orders m- 



n> 1, whereas ctoo, 0^10 and croi consist of parts belonging 
to (T and (T. The expansion of f comprises the three terms 
foo, fio and foi according to ( |38| ) and (p6|), and therefore 
fmo in ( |63|) is a completely different fictitious body force. 
For these reasons, the "-notation of the projection method 
is not appropriate in case of Eq. (|6^). 

As an important result of the above considerations, 
the perturbation expansions of the elastic fields can be 
calculated by solving the equations of balance with the 
differential operator Doo(V) instead of the full D(V), a 
task which proves to be relatively easy. We have per- 
formed our calculations with the help of M aple V |^, and 
made use of the method described in Section |3.3| . There- 
fore, all displacement fields are compounded by two parts: 
7m« = 7mn + 7'm„ for all mn. 



3.5 Solving the recursion formulae 

The calculation of the displacement fields 7^ 
use of the fact that they depend on only: 



7™„(x) = 7mn(<^) : 



can make 



(64) 



with each component 7^ 



qWP being a real 



combination of complex harmonics. Therefore, things are 
much easier when expressing all derivatives ( pq ) in terms 
of d/dr, d/dij) and d/dz, with d/dr and d/oz vanishing 
when applied to the displacement fields (pi]). 

^ MapleV Release 4 ©1981-1996 by Waterloo Maple Inc. and 
MapleV Release 5.1 ©1981-1998 by Waterloo Maple Inc.. 
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Transforming Doo(V) yields the equations of balance 
of zeroth order in the form 



-(A + 3m)^ 
-Ai+(A + 2/i)^ 











+ 





-0, 



0, 



M5i< + /r = 0, 



(65) 



where k £ z}. From (|45|), all possible (div . . . )\x^o 

in Eqs. (61), ( |63| ) and in the analogous equations for other 
orders depend on r like r~^, and in ( |65| ) we use the nota- 
tion (div . . . )|x^o — -pii, where the length force f depends 
on (j) only. Eqs. ( |65| ) justify the ansatz (|6^). The order in- 
dices are omitted here. 

In the second step of the calculation of 7„„, the diver- 
gence of the stress tensor on the dislocation line must be 
brought to zero. This leads to the additional displacement 
field entering 7„„. To perform this second step, the 
fields 7^ resulting from (^ with the prototypical forces 
i'l^ defined by f'l^ ^ = Sap must be available. We have 
calculated these 7^ with the help of ( |2^ ) for zeroth order. 
They are given in Appendix B. From there it is obvious 
that Inr-terms and hereby dependencies on r join the dis- 
placement fields. 

For the purpose of this paper, we have calculated the 
elastic displacement, strain and stress fields of two- and 
threefold dislocations up to orders mn where m + n < 5, 
and for fivefold dislocations up to m + n < 15. 

For fivefold dislocations and vanishing phonon-phason- 
coupling, the elastic fields are known in closed form p6t . 
When we neglect the phonon-phason-coupling in our al- 
gorithm, then our analytical expressions agree with an ex- 
pansion of this closed solution into a perturbation series in 
/X4 — 115 up to all calculated orders. But also our algorithm 
clearly shows, that neither a closed solution is possible for 
finite phonon-phason-coupling nor for other directions of 
the dislocation line. 

The solution of [26] has been obtained by the more 
frequently used generalized Green's method. The calcu- 
lated displacements are applied to atomic models of dis- 
locations along a fivefold direction of i-AlPdMn to obtain 
dislocat ion- induced atomic positions. As we shall prove in 
Section |4.2| , the true fivefold symmetry of the displace- 
ment fields is exclusively induced by the phonon-phason- 
coupling, but investigating the effect of the phonon-phason- 
coupling on the atomic positions is beyond the scope of 
this paper. 



4 Results and discussion 
4.1 General results 

The strain and stress fields of all orders behave like - , with 
r being the distance from the dislocation line. This is ex- 



actly compatible with what one would expect from the ex- 
act solution, provided e.g. by the projection method, and 
is also compatible with results from numerical methods. 

Apart from zeroth order, the displacement fields show 
positive parity perpendicular to the dislocation line: 



(66) 



In contrast to this, the strain and stress fields of all orders 
exhibit negative parity: 

^mn( ^) — ^m7i(^) 5 ^rnni, ^) — ^miiip^^ • i^^) 

These properties follow from the equations of balance (|l^ ) 
and (p9[), respectively, when using the fact that all ficti- 
tious forces occuring in the projection method have posi- 
tive parity. Note that (|6^) guarantees for the total torque 
acting on the dislocation core to vanish. 

The components em„^33 and em„.33 of the strain ten- 
sors vanish. This is a consequence of the infinite geom- 
etry. Except for the case of twofold screw dislocations, 



"mn, 33 



^mn, 33 = is uot Compatible with cr^ 



mn,33 





and (T^„ 33 = 0. Therefore, stresses are present parallel to 
the dislocation lines. It is not quite clear how to interpret 
these statements in phason space. 

The simplest dislocations are the twofold ones. Twofold 
edge dislocations induce only displacements perpendicular 
to the dislocation line, and screw dislocations only parallel 
to the dislocation line. In case of the other line directions. 
Burgers vectors and displacements perpendicular and par- 
allel to the dislocation lines are no longer decoupled. 

For every order mn where m is even, the phononic 
displacement, strain and stress fields depend linearly on 
b" only, whereas the phasonic fields depend on only. 
The contrary is true for any order mn, where m is odd. 
Additionally, we have Uo„ = 0, Sq^ = and cTq^ = for 
all n > 0. 

In icosahedral quasicrystals, it is, in general, not pos- 
sible to compute the elastic fields belonging to a Burgers 
vector arising from b by rotation about the dislocation 
line from the fields belonging to b by the same rotation. 
This method can, for example, be applied to pentagonal 
quasicrystals 0]. 



4.2 Displacement fields of fivefold dislocations 

The following performanc es r efer to the fivefold coordi- 
nate systems introduced in 3.2. The classification into pha- 



sonic edge and screw dislocations is problematic in some 
respects [^|. To our opinion, this discrimination makes 
sense, because only phasonic Burgers vectors which are 
parallel to the z-axes of the coordinate systems of 3^ do 
not break the n-fold symmetries induced by phononic n- 
fold screw dislocations. 

The higher displacement corrections shown in the next 
two subsections are computed on the circle i/a^^ + y^ = 
10 b in the xy-plane. Here, b is the length of the different 
Burgers vectors we have assumed (see Eqs. (|8|), (|69|)). 
We have chosen this circle, because for our choice of elas- 
tic constants, this is the lower limit for no strain tensor 



M. Ricker, J. Bachteler, H.-R. Trebin: Elastic theory of icosahedral quasicrystals - application to straight dislocations 9 



component to exceed an absolute value of 0.1, which we 
regard as limit for the validity of linear elasticity. In the 
diagrams with displacement components Uz parallel to the 
dislocation line, the label L stands for the arc length of 
this circle. All displacement fields are standardized to have 
a vanishing net displacement on the circle, which has been 
achieved by adding suitable constant displacement fields. 

For the calculation of the elastic fields, the elastic con- 
stants of i-AlPdMn have been used. The phononic elas- 
tic constants are A = 85GPa and fj, = 65GPa For 
the phasonic elastic constants, we have taken the val- 
ues Ki/kT = 0.1/atom and K2/kT = -0.05/atom and 
supposed a quenching temperature of 500° C This 
leads to our elastic constants ^4 — 0.012 GPa and ^5 = 
0.12 GPa. The phonon-phason-coupling has been supposed 
to be /L(3 = 1 GPa. 

Note that the displacements alone are no good quanti- 
ties to measure the elastic deformation qualitatively, but 
the strains are. Note also, that the phononic displacement 
corrections are very small compared to the phononic Burg- 
ers vectors, and that all the corrections shown here vanish 
with the phonon-phason-coupling /13. 



4.2.1 Fivefold edge dislocation 

We have chosen a six-dimensional, prototypical Burgers 
vector 

b"'"" = [0, 6, 0]' . (68) 

For the length of the phononic part b" , 1 A is the right or- 
der of magnitude, whereas under plastic deformation the 
phasonic part b™ increases and becomes up to 100 times 
larger than b" Q . Pure edge or screw dislocations appear 
very rarely in reality. Due to the edge type of the disloca- 
tion defined by (^8|), the fivefold symmetry of the disloca- 
tion line is broken, as can be seen from the displacement, 
strain and stress fields. 

Figs. ^ and ^show the xy-components of the phononic 
displacements of order mn = 10 and mn — 20. They con- 
tribute most to the total correction of the phononic xy- 
displacement field, which is shown in Fig. |5|. Obviously, 
the superposition of these two low orders already provides 
a good approximation for the full correction. 

In Figs. H and 0, the z-components of order mn = 20 
and of all orders mn with m + n < 15 are shown. Again, 
the orders mn — 10 and mn = 20 give the essential 
contribution to the full displacement field. In zeroth or- 
der, no z-component is present. The z-component of order 
mn = 10 is not shown here. It proves to have a negative 
cos 20-dependence, which in 2000-fold magnification has 
a similar absolute amplitude as the components of Fig. ^. 

The perfect twofold symmetry of the z-components is 
remarkable. This is due to the condition (|66|), which the 
displacement fields have to fulfill. Note that the compo- 
nents of Figs. fulfill ( |66| ) also. In addition to this sym- 
metry, the displacement fields show mirror symmetry with 
respect to the xz-plane. As an evident explanation, the 




Fig. 3. Phononic a;y-displacement field of the edge disloca- 
tion (^), order mn = 10. Magnification of the displacement 
vectors; 4000. 



y + 




Fig. 4. Phononic xy-displacement field of the edge disloca- 
tion (p^), order mn — 20. Magnification of the displacement 
vectors: 6000. 



Volterra cut procedure ||T^,Q should be performed sym- 
metrically with respect to a cut in the physical a;z-plane, 
so that this plane remains a mirror plane even after the 
introduction of the dislocation (see Figs. |l|, ||). Our results 
show that the phasonic fields, which have to be subjected 
to the associated orthogonal space symmetry operation, 
also have this symmetry, as it must be. 

The volume contraction = tr (e") and the eigen- 
values of the phononic stress tensor cr" of zeroth order do 
not change substantially from the contribution of higher 
orders. 

When considering phasonic components, the typical 
displacements are more than a hundred times larger than 
their phononic counterparts. The typical phasonic strains 
result to be ten times larger than the phononic strains. 
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y f 




Fig . 5. Phononic xy-displacement field of the edge dislocation 
(pq), all orders mn with m + n < 15 except zeroth order. 
Magnification of the displacement vectors: 2500. 




Fig . 8. Phononic a;j/-displacement field of the screw dislocation 
( |69[ ), order mn — 20. Magnification of displacement vectors: 
2500. 








lOTtb 



Fis 



6. Phononic ^-displacement field of the edge dislocation 
order mn = 20. Magnification of the arrows: 1200. 




lOjib 



Fig . 7. Phononic z-displacement field of the edge dislocation 
(pSj), all orders mn with m + n < 15. Magnification of the 
arrows: 2000. 



We have investigated the convergence behaviour of our 
perturbation expansion and also calculated the exact dis- 
placement fields by the numerical Eshelby's method 
2l|. From there, the correction displacement fields of Figs, 
^nd are very close to the exact ones. The same is true 
for the correction displacement field for the fivefold screw 
dislocation, which is discussed below. 




Fig. 9. Phononic a;j/-displacement field of the screw disloca- 
tion ([69[), all orders mn with m + n < 15. Magnification of 
displacement vectors: 1500. 



4.2.2 Fivefold screw dislocation 

The Burgers vector of a screw dislocation has been chosen 
to be 



b"'™ = [0,0,5]* 



(69) 



A screw dislocation induces no symmetry breaking paral- 
lel to the dislocation line. Therefore, the fivefold symme- 
try should be recognizable from the phononic as from the 
phasonic displacement fields. 

From Fig. the fivefold symmetry of the phonon dis- 
placement field is obvious. Again, mn = 20 is the or- 
der which provides the largest contribution to the com- 
plete xy-displacement field, which is shown in Fig. |^. Al- 
though the displacement fields look qualitatively very sim- 
ilar, slight differences are present in their functional de- 
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5b 



^ -5b 



-5 b 



1:. 



5b 



Fig. 10. Phasonic a;y-displacement field of the screw disloca- 
tion (ji^), order mn — 10. Magnification of the displacement 
vectors: 6. 



mn = 40, but the z-componciit of order mn = 40 is about 
11 times as large as in order mn = 31, and provides the 
largest z-correction. The graph of Fig. In] is not exactly 
symmetric with respect to its extremal points. 

The displacement fields exhibit no invariance under 
the action of an appropriate mirror plane containing the 
z-axis, as in case of the edge dislocation (^). 

The eigenvalues of the stress tensor cr" , which are con- 
stant in zeroth order, become clearly fivefold symmetric 
when taking into account higher orders. The main change 
away from isotropic behaviour is provided by orders mn 
where m + n < 2. Higher orders provide also a fivefold 
volume contraction, whereas in zeroth order there is none. 
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Fig . 11. Phononic ^-displacement field of the screw dislocation 
(p9|), all orders mn with m + n < 15 except zeroth order. 
Magnification of the graph: 400000. 



Appendix A: Icosahedral irreducible strains 



pendence on x. No z-component is present in order mn = 
20. 

In case of phason displacements, the fivefold symme- 
try is not easy to recognize. Fig. |l^ shows the phasonic 
displacement field of order mn = 10, which provides the 
largest phasonic xy-displacements occuring in any order. 
Again, no z-component is present. Note the small magni- 
fication, indicating the magnitude of the phason displace- 
ments, which are appreciable compared with the length 
of the phasonic Burgers vector. The position vectors used 
in Fig. ^ refer to the physical coordinate system of Fig. 
H, but the components of the displacement vectors belong 
to the phason coordinate system of Fig. ^. Looking at the 
indices of the projected hyperspace basis vectors in 
Figs. |l| and |[ it is clear, that a rotation about angle ^ 
in physical space corresponds to a — ^-rotation in phason 
space. Therefore, fivefold invariance of the displacements 
in Fig. |l^ refers to a —^-rotation of the vectors and a 
^-rotation of the positions. 

The total correction of the displacement field has a 
z-component which is shown in Fig. |ll|. Its tenfold sym- 
metry is surprising, but it is the only way for the system 
to obey the constraint (jo^) . The lowest orders having such 
a tenfold z-component in phonon space are mn = 31 and 



The icosahedral irreducible strains given here are taken 
from m. 



1 

7! 



£a = 



+£33) 



2^3 



11 



+ 7^e22 + (^+ 7)^33)' 



- ^^22 

\/2e5'2 



+ £33) 



-23 
.u 
31 



1 

71 



1 

7i 



22 



'33 



32 



- Tel 

- TS 



23 



- re 



31 



- 2e 



(70) 



r^lS 

V3(eT'i 
^u+e^2 

V2(re^2-7£^3 



33 
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Appendix B: Elastic displacement fields of ze- 
roth order 

From (Q), (^Tj) and (|62|), the displacement field of zeroth 
order is 



_ 1 



+ b. 



r ^ (A + 2^)(a;2+y2)/ 



X + 2fj. 7-0 (A + 2^)(x2 + ?;2). 
(A + fi)xy \ 



(A + 2^)(a;2+y2)J 



X + 2fj. 7-0 (A + 2^)(x2 + ?;2) 



5", 



(71) 



Here, 700,1 = uqq.x, ■ ■ ■ ,700,6 = woo,z- This displacement 
field is identical to the solution presented in when 
K2 — 0. K2 = corresponds to identical elastic constants 
fi4 — /i5, what is obvious from ([74[). The factor rg is in- 
troduced for dimensional reasons, and sometimes it is in- 
terpreted as inner cut-off for the linear elasticity Q| . 

The application of the projection method requires the 
knowledge of the solutions of the elastic equations of bal- 
ance 



Doo(V)7::+'5(x)fL-0, 



(72) 



where the forces 6{x)i'l^, f'l^p = 6ai3, a, = 1, . . . , 6, 
represent point-like stress sources on the dislocation line 
(these prototypical forces are not body forces!). We have 
obtained the resulting fields using Eq. (23) with Goo(x — 
x'). The 7^ are 

(A -I- /i) a;2 



7'/ 



I2 



13 



74 



1'^ 



76 



1 

An 



Ai(A + 2/i) ro fi{X + 2fi){x^ +y^)' 



(A + n) xy 



1 

47r 



1 

2^ 
1 

2^ 
1 

2^ 
1 

2^ 



MA + 2M)(x2+t/2) 
(A + ^) xy 



0, 0, 0, 



m(A + 2/^)(x2 -l-?;2) ' ^(A + 2^) 



J^o 



(A + A^)?;^ 



/Lt(A-h2/i)(a:;2 -t-y2 

0, 0, -iln — , 0, 0, 

1 r 

0, 0, 0, In — , 0, 

fo 



, 0, 0, 0, 



0, 0, 0, 0, - 
0, 0, 0, 0, 



1 r 
— In — , 

fo 



(73) 

These fields don't have the unit of a length, but our algo- 
rithm multiplies them by appropriate factors to become 
displacement fields. 



Appendix C: Elastic constants 

Different notations concerning the five independent elastic 
constants of icosahedral quasicrystals are in use. 

To every elastic constant belongs a quadratic invariant 
of strain components. Different sets of elastic constants 
can be compared with each other, when the respective 
invariants are formulated in one and the same coordi- 
nate system. This can be achieved by means of coordi- 
nate transformations. Relations between different elastic 
constants arise from expressing one set of invariants by 
another and comparing the coefficients. 

All authors we adress also use the two Lame-constants 
A and /i (|l^), but a different phonon-phason-coupling and 
other phasonic elastic constants. 

In PI and [ p6| , the phonon-phason-coupling is de- 
noted i?, and the two phasonic elastic constants are Ki 
and K2. The coordinate systems in phonon and phason 
spaces in these papers are transformed into our fivefold 



if. 



when substituting x —x, y —y, z 



and vice 



versa. Transforming our elastic tensor C5 in this manner 
allows a comparison with the parts R and K of the elastic 
tensor given explicitly in ||l^. The relations between the 
elastic constants are 



fJ-3 = \/6i?, 

^4 - if 1 - 2K2 ., 

fi5^Ki+K2, 



Ki = -(m4 + 2^15) , 

K2 = --^(^4 - Ms) 



(74) 



The phonon-phason-coupling of |2^ is denoted K3 and 
the phasonic elastic constants are Ki and if 2 , but with an- 
other meaning as in [HI and ||2^. To compare the elastic 
energy expression of [ 27] with our Eq. (|lj) , one must sub- 
stitute £^ by for all This is because the definition 

^ij ^ lix^ [13 different from our definition (|8|). In a 
second step, this new expression must be transformed into 
our coordinate systems of Figs. |l|, ^ what is achieved by 
subjecting all components e^^^ to the coordinate trans- 
formation X y, y ^ —X, z ^ z. From the elastic energy 
density in this new form, we see that the elastic constants 
are related by 



V6if3, 



1 



V6 



Ma : 



At4 = if 1 + gif2 

4 

= ifi - 2if2 



ifi = -(4/i4 + 5M5), (75) 
if2 = -^(^4 - Ais) • 
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